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Abstract 

We derive the nonequilibrium real-time evolution of an 0{N) — 
invariant scalar quantum field theory in the presence of a nonvanishing 
expectation value of the quantum field. Using a systematic 
expansion of the 2PI effective action to next-to-leading order, we 
obtain nonperturbative evolution equations which include scattering 
and memory effects. The equivalence of the direct method, which 
requires the resummation of an infinite number of skeleton diagrams, 
with the auxiliary-field formalism, which involves only one diagram at 
next-to-leading order, is shown. 
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1 Introduction 



Nonequilibrium quantum field theory has a wide range of applications, 
including current and upcoming relativistic heavy-ion collision experiments 
at RHIC and LHC, phase transitions in the early universe or the formation of 
Bose-Einstein condensates in the laboratory. Important theoretical progress 
has been achieved with effective descriptions based on a separation of scales 
in the weak-coupling limit [Q, or for systems close to equilibrium using 
approximations such as (non)linear response or gradient expansions 0. 
However, the description of far-from-equilibrium dynamics is still in its 
infancies. The situation is complicated by the fact that typically there is 
no clear separation of scales which is valid at all times and it is often difficult 
to identify a small expansion parameter. For example, standard perturbation 
theory is plagued by the problem that a secular (unbounded) time evolution 
prevents the description of the late-time behavior of quantum fields. 

Practicable nonperturbative approximation schemes may be based on the 
two-particle irreducible (2P/) generating functional for Green's functions @ . 
Recently, a systematic expansion of the 2PI effective action has been 
proposed and applied to a scalar 0(A^)-symmetric quantum field theory Q. 
The approach provides a controlled nonperturbative description of far-from- 
equilibrium dynamics at early times as well as the late-time approach to 
thermal equilibrium and can be applied in extreme nonequilibrium situations 

0. This is in contrast to the standard expansion of the IP I effective 
action which is secular in time once direct scattering is taken into account 
0. The 2PI-1/N expansion extends previous successful descriptions of the 
large-time behavior of quantum fields 0, §], which employ the loop expansion 
of the 2PI effective action relevant at weak couplings ^, |TU|. For systems 
in or close to equilibrium recent applications of the loop expansion can be 
found in Refs. [0, [T^, |TB|, n|. A expansion has the advantage over the 



loop expansion that it is not restricted to small couplings, an observation 
made in the context of nonequilibrium quantum field dynamics already for 



quite some time However, in order to describe quantum scattering and 
thermalization the inclusion of the usually discarded next-to-leading order 
(NLO) contributions is crucial. 

In Ref. [0 the 2PI-1/N expansion has been carried out to NLO in the 
symmetric regime for a vanishing expectation value of the quantum field. 
Here, we derive the evolution equations from the NLO approximation of the 
2PI effective action for nonzero 0. A nonvanishing field expectation value is 
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important to describe the physics of heavy ion colhsions, where the presence 
of a substantial scalar quark-antiquark condensate signals the spontaneous 
breakdown of chiral symmetry. (In the language of the 0(4) linear sigma 
model for two-flavor QCD one has ~ {qq)-) Important nonequilibrium 
applications in this context include the formation of disoriented chiral 



condensates (DCCs, cf. [T^) or the decay JT^ of parity odd metastable states 
in hot QCD [T^. Similar far-from-equilibrium applications in inflationary 
cosmology concern the phenomenon of preheating, where the dynamics of 
the infiaton is expected to trigger explosive particle production (see, e.g., 



T9|). The availability of a quantum field description of the dynamics is 
important for cases where the effectively classical approximation (in the 
context of inflation, see e.g. PU[) i^^Y iiot be appropriate. 



In this work, we perform a systematic expansion of the 2PI effective 
action to next-to-leading order. A detailed discussion of the classification of 
diagrams is given in Sec. ^ This allows us to give the effective action (Sec. ^ 
and the equations of motion (Sec. ^ in a straightforward manner. In Sec. ^ 
it is shown how this result is obtained using the auxiliary-field formalism 



21| . The causal equations for the spectral and statistical functions are listed 
in Sec. |^ and a perturbative approximation to second order in the coupling 
constant is given in Sec. H. In two Appendices we discuss further features of 
the equations and the realization of Goldstone's theorem. 



2 The 2PI effective action 

We consider a real scalar A^-component quantum field (fa (a = 1, • • • , A^) 
with a classical 0(A^)-invariant action 



X 



(2.1) 

Summation over repeated indices is implied and x = (x°,x). All correlation 
functions of the quantum theory can be obtained from the effective action 
r[0,G], the two-particle irreducible (2P/) generating functional for Green's 
functions parametrized by the macroscopic field (pa{x) and the composite 
field Gab{x,y), given by 

0,(x) = (<^„(x)), (2.2) 
Gab{x,y) = {Tcipa{x)(pb{y)) - {Mx)){(pb{y)). (2.3) 
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The brackets denote the expectation value with respect to the density matrix 
and Tq denotes time-ordering along a contour C in the complex time plane. 
At this stage the explicit form of the contour is not needed. 

A discussion of the defining functional integral of the 2PI effective action 
can be found in Ref. f^. Following it is convenient to parametrize the 
2PI effective action as 

r[0, G] = S[4>] + '-Tt InG-' + ^Tr Gq\^) G + Tsl^, G] + const. (2.4) 
Here the classical inverse propagator iGQ^{(f)) is given by 



6(f)aix)6(f)biy) 



Sab + ] ^cix - y) (2.5) 



with 0^ = (pa(pa and 6c{x — y) = 6c{x^ — y°)5'^(x — y). The contribution 
r2 [4>, G] is given by all closed 2PI graphs[| with the propagator lines set 
equal to G 0]. The effective interaction vertices are obtained from the terms 
cubic and higher in (p in the classical action (|2.1| ) after shifting the field 
(j) + (p. In presence of the nonvanishing expectation value (p this results 
in a cubic and a quartic vertex 

Cintix; 0, if) = —^(l)aix)ipaix)ipi,ix)ipbix) - ^ (x)(/?a (x)]^ . (2.6) 

Dynamical equations for (pa and Gab can be found by minimizing the 
effective action. In the absence of external sources physical solutions require 

d(pa{x) 

which leads to the macroscopic field evolution equation 

+ + [(l)'^{x) + Gbb{x,x)]^ (paix) = -^(pb{x)Gba{x,x) 



S(pa{x) 



graph is two-particle irreducible if it does not become disconnected upon opening 
two propagator lines. 
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as well as 



which leads to 



with 



6Gabix,y) 
j:ab{x,y;(f),G) = 2%- 



(2.9) 



(2.10) 



(2.11) 



6Gab{x,y)' 

Eq. (|2.10| ) can be rewritten as a partial differential equation suitable for 
initial value problems by convolution with G, 



<^o.afe(^> z)Gbc{z, y) = ^ac{x, z)Gcb{z, y) + 6ab5c{x - y) , (2.12) 



where the shorthand notation = dz^ J dz is employed. With the 
classical inverse propagator (|2.5|) this differential equation reads explicitly 



Gabix,y) = ^(j)aix)(f)c{x)Gcbix,y) 
+ i I Eac{x,z;(f),G)Gcb{z,y) +i6abSc{x -y). (2.13) 

J z 

The evolution of (pa and Gab is determined by Eqs. ( p.8| ) and ( p.l3| ), once 
r2 [0, G] is specified. 



3 The 2PI-1/N expansion 

In this section we discuss the 2PJ-1/A^ expansion proposed in Ref. in more 
detail. We present a simple classification scheme based on 0(A^)-invariants 
which parametrize the 2PI diagrams contributing to r[0, G]. The effective 
action at a given order in 1/A^ can then be obtained from a straightforward 
summation of the diagrams contributing at that order. The NLO result 
for r[</), G] for general field configurations (0, G) is presented in Sec. ^. 
Below we will compare the direct summation procedure described here with 
an alternative method, which employs an auxiliary field to simplify the 
summation (cf. Sec. H). 
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3.1 Counting rules for irreducible 0(A^)-invariants 



Following standard procedures ||22| the interaction term of the classical action 
in Eq. (2J.) is written such that S[(j)] scales proportional to N. From the 



fields (pa alone one can construct only one independent invariant under 0{N) 
rotations, which can be taken as tr (jxp = 0^ = (j)a4'a ~ N. The minimum 0o 
of the classical effective potential for this theory is given by 0q = N{—Qm?/ A) 
for negative mass-squared and scales proportional to N . Similarly, the 
trace with respect to the field indices of the classical propagator Go is of 
order A^. 

The 2PI effective action is a singlet under 0{N) rotations and 
parametrized by the two fields (pa and Gab- To write down the possible 
0{N) invariants, which can be constructed from these fields, we note that 
the number of 0-fields has to be even in order to construct an 0(A^)-singlet. 
For a compact notation we use {(pcfyab = 4>a4>b- Consider a general singlet 
composed of arbitrary powers of {(p(p)ab and Gab (we neglect all space-time 
dependencies) 



s = tr[(#)P^G«V..(#)P"G'?"] 

= (02)Pi+-+P"-"tr(00G'^i) ...tr(00G«") (3.1) 

for Pi 7^ 0, where the pi and qi {i = 1, . . . ,n) are any positive integer. (For 
all Pi = the RHS is given by ti G'^^^ "^'^" .) The second line in the above 
equation follows from simple contraction of the field indices and corresponds 
to the fact that all functions of and G, which are singlets under 0{N), 
can be built from the irreducible (i.e. nonfactorizable in field-index space) 
invariants 

02, tr(G") and tr(00G'"). (3.2) 

We note that for given N only the invariants with n < N are irreducible. 
However, we will see below that for lower orders in and for sufficiently 
large one has n < N. In particular, for the next-to-leading order 
approximation one finds that only invariants with n < 2 appear. 

Since each single graph contributing to r[0,G] is an 0(A^)-singlet, we 
can express them with the help of the set of invariants in Eq. ( |3.2|) . The 
factors of in a given graph have two origins: each irreducible invariant is 
taken to scale proportional to since it contains exactly one trace over the 
field indices, while each vertex provides a factor of 1/A^. The leading order 
(LO) graphs then scale proportional to A^ (as the classical action S), next- 
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Figure 1: Graphical representation of the 0-dependent contributions for 
r2 = 0. The crosses denote field insertions ~ (pafpa for the left figure, which 
contributes at leading order, and ~ 4>a<Pb for the right figure contributing at 
next-to-leading order. 

to-leading order (NLO) contributions are of order one, the next-to-next-to- 
leading order (NNLO) scales as 1/iV and so on. This provides a well-defined 
expansion of r[0, G] in powers of 1/N. We stress that by construction each 
order in the 2P/-1/A^ expansion respects 0{N) symmetry. In particular, this 
is crucial for the validity of Goldstone's theorem in the case of spontaneous 
symmetry breaking. For constant 7^ one observes that the mass matrix 
~ (52r[0,G'('^*^*)(0)]/50a506, with G(^*'^*)(0) being the solution of Eq. 
contains (A^— 1) massless "transverse" modes and one massive "longitudinal" 
mode (see Appendix A for a discussion). 

3.2 Classification of diagrams 

In the following we will classify the various diagrams contributing to r[0, G]. 
The expression ( |2.4| ) for the 2PI effective action contains, besides the classical 
action, the one-loop contribution proportional to Tr lnG~^ + TrGQ^(0)G 
and a nonvanishing r2[0, G] if higher loops are taken into account. The one- 
loop term contains both LO and NLO contributions. The logarithmic term 
corresponds, in absence of other terms, simply to the free field effective action 
and scales proportional to the number of field components A^. To separate 
the LO and NLO contributions at the one-loop level consider the second term 
TrGQ^(0)G. From the form of the classical propagator ( |2.5|) one observes 
that it can be decomposed into a term proportional to tr(G) ~ A^ and terms 
(neglecting all the space-time structure) 

~74N#)tr(G)+2tr(#G)]. (3.3) 
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a b 

oo 

a b 



Figure 2: LO contribution to the 2PI effective action. 



This can be seen as the sum of two "2P/ one-loop graphs" with field insertion 
~ (j)a4>a and ~ 0a0b, respectively. Counting the factors of N coming from the 
traces and the prefactor, one observes that only the first of the two terms in 



( |3.3| ) contributes at LO, while the second one is of order one. 

We now turn to r2 [(f), G] which contains all 2PI diagrams beyond the 
one-loop level. The graphs are constructed from the three-point vertex and 
the four-point vertex in the interaction Lagrangian (|2.6|) with G associated 
to the propagator lines. From the three- vertex ~ 0^ one easily observes that 
one cannot construct a diagram which has a field insertion ~ 0^ = (f)a4'a- As 
a consequence, all loop diagrams beyond the one-loop level can only depend 
on the basic invariants tr(G"') and tr((/)0G'"). Furthermore, we note that 
the invariant tr(G) = Gaa can only appear if the graph contains a tadpole as 
shown in Figs. ^ and|^. The first graph. Fig. |^, is proportional to the product 
of the term [tr(G')]^ and a factor from the quartic vertex. It therefore 
contributes at LO to the effective action. The graph in Fig. |^ is proportional 
to ii{G'^)/N and of order one. The two-loop graphs shown in these figures 
are indeed the only two-particle irreducible graphs which contain a tadpole. 
As a consequence, the only invariants that can arise beyond two loops are 
tr(G"^2) and tr(00G"^i). 

Using these considerations, one can now straightforwardly continue to 




Figure 3: NLO "double bubble" contribution. 
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Figure 4: NLO 0-independent contribution to the 2PI effective action. 
Higher loop diagrams in the infinite series can be obtained from the previous 
one by introducing another "rung" with two propagator hues at each vertex. 
The resummed series including the prefactors not displayed in the figure is 
given by the first term in Eq. ( |4.3| ). 



classify all the diagrams at a given order of the expansion. Consider a 
graph with V3 three-point and V4 four-point vertices. The number of internal 
lines / is determined by the familiar relation 



21 = 3V3 + 4^4- 



(3.4) 



We observe again that V3 has to be even. This graph contains V3/2 field 
insertions ~ 00, 2V4-I-IV3 propagators G and comes with a factor {X/N)^^~^^* 
from the vertices. The highest power of is obtained by contracting the 
and the G field such that the largest number of invariants tr(G"'^^) and 
tr(00G'"'^^) appears. The structure with the highest power of iV, 



A 

N. 



V3+V4 



tr(00G) tr(G 



V4+ 



V3 



(3.5) 



is of order one and contributes therefore at NLO. 

For V3 = the corresponding diagrams are shown in Fig. § [Q. For 
V3 = 2 the infinite series of graphs is presented in Fig. ^. We will now 
argue that graphs with V3 > 4 are not two-particle irreducible such that we 
have classified the complete NLO contribution. We note that the invariant 
tr(00G) connects the field insertion ~ 0a06 with a single propagator line. 
If a diagram contains that invariant more than once it is always possible to 
disconnect the graph by cutting two lines.0 Consequently, a diagram with 
the above structure and with V3 ^ 4 is not two-particle irreducible and does 



^ To see this, note that the number of loops L in these diagrams is given by the standard 
relation 



L = I-Vs-V4 + l^V4 + 



1, 



(3.6) 
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Figure 5: NLO 0-dependent contribution to the 2PI effective action. Each 
diagram in the infinite series can be obtained from the previous one by 
introducing another "rung" with two propagator hues at each vertex. The 
resummed series is given by the second term in Eq. ( [4.3|) . The complete NLO 
contribution is given by the sum of the diagrams in Figs. ^ and |^. 



not contribute to F [(/),(?]. We therefore have classified all possible diagrams 
which contribute to the 2PI effective action at NLO and present the result in 
the next section. In Appendix A we discuss possible further approximations 
consistent with an expansion in powers of 



4 The 2PI effective action at NLO 



We write 

F2[0,G]=F^«[G]+Ff«[0,G] + ... (4.1) 

where F2° denotes the leading order (LO) and F^'"^ the next-to-leading order 
(NLO) contributions. The LO contribution to ^2[G] is given by the diagram 
presented in Fig. ^ 

r^°[G'] = I^Gaa{x,x)G,,{x,x). (4.2) 

This contribution is 0-independent. 

The NLO contribution consists of an infinite series of diagrams which fall 
into two classes. The first class is independent of (p and is constructed with 
only quartic vertices. It contains the complete NLO contribution when = 
and has been resummed in Ref. |Q. The diagrams are shown in Fig. ^ The 



where we have used Eq. (3.4) for the second equaUty. Tadpole diagrams are forbidden in 
those graphs and each "bubble" tr(G'^) corresponds to one closed loop. From the above 
relation one then observes that the total number of loops in the diagram is given by the 
number of bubbles plus one. In particular, one can disconnect the diagram by cutting two 
lines connecting the fields insertions (the terms ^ tr(0(/)G)). 
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second class depends on and contributes for the case of a nonzero field 
expectation value. This series of diagrams is shown in Fig. |] and can be 
resummed as well. As a result, the NLO contribution to T2[4>, G] is given by 
the sum of the two resummed expressions and readsQ 

rf''[0,G] = ^Tr,Ln[B(G)] + ^^ I{x,y;G)U^)G,b{x,y)My)- (4-3) 

In the above equation we have defined 

B(x, y; G) = 5cix - y) + i-^Gabix, y)Gabix, y), (4.4) 



and the logarithm in Eq. ( [4.3|) sums the infinite series shown in Fig. 
Tr^Ln[B(G')] = j (^i-^Gab{x,x)Gab{x,x) 

~\ / (^6^*^'^**^^'^)^'^''*^^'^)) (^^-^^<^'b'{y^^)Ga'b'{y,x) 
+ ... (4.5) 
The function I(a;, y; G) is defined by 

A A /" 

I(x, y\ G) = —Gabix, y)Gab{x, y) - i— / I(x, z; G)Gab{z, y)Gab{z, y) , 

(4.6) 

and resums the infinite chain of "bubble" graphs, which can be seen by re- 
expanding the series. The functions I(a;, G) and the inverse of B(x, G) 
are closely related by 

B-i(x, y- G) = 5c{x - y) - il{x, y; G) , (4.7) 

which follows from convoluting Eq. ( [4.4|) with B~^ and using Eq. ( [4.6|) . We 
note that B and I do not depend on 0, and T2[4>, G] is only quadratic in at 
NLO. Hence, the complete effective action at NLO contains only quadratic 
and quartic terms in 0. 

^ Besides the dynamical field degrees of freedom (j> and G we will introduce a number of 
quantities which are (resummed) functions of these fields. These functions will be denoted 
by either boldface or Greek letters in the following. 
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5 The equations of motion 

From the 2PI effective action r[0, G] at NLO we find equations of motion 
for tlie macroscopic fields and G, as indicated in Sec. ^. Tlie equation for 
tlie field expectation value ( |2.8| ) reads at NLO 

-(^^ + w? + -^ + Gcc{x, x)] ^ = Ka(x, x). (5.1) 

We have written the LO contribution of the evolution equation on the LHS 
and combined the NLO contribution a^ 

K,(x, y) ^ K„(x, y;^,G)^^J^ B-\x, z; G)Gabiy, z)M^), (5.2) 

evaluated at x = y. We have written Ka(x, ?/) as a function of x and y for 
later convenience. 

Eq. (|2.10|) for G~f^ (x, y) can be completed using Eq. (|2.11|) for the self- 
energy. The LO contribution is simply 

^^hi.^^y) = -i-^Gccix,x)6ab5cix - y). (5.3) 

To obtain the NLO contribution the following identity may be helpful 

Sl{u,v;G) .6B-\u,v;G) A , , 

aTw V = ^ xr^ ( ^ = {u,x;G)Gab{x,y)B {y,v;G), 

dGab{x,y) dGab{x,y) 3A/ 

(5.4) 

where we used that 

6B{u,v;G) • A , . , . , . , . 

-TTT-f r = ^:T^Gab{x, y) dc{u - x)dc{v - y). (5.5) 

bGab[x,y) SN 

Collecting all the pieces, we find that Eq. ( p.lOj ) for the inverse propagator 
can be written as 

iG~,^{x, y) = - {Ux + + ^ W{x) + Gcc{x, x)] ^ 5abbc{x - y) 

^ ^'^{x,y]G)(t)a{x)(l)b{y) + iD{x,y)Gab{x,y), (5.6) 



^Note that the classical inverse propagator iGq ^ contains a LO and a NLO part. 
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with the definition 

D(x, y) = D(x, y- 0, G) = G) 

+ {^) I B-\x,u-G)Mu)Gab{u,v)Mv)B'\v,y-G). {5.7) 

Since is of order one and D of order the first fine on the RHS of Eq. 
(|5.6| ) corresponds to the LO and the second fine to the NLO contribution. 

Eqs. (13, |5]|) together with Eqs. (|3, ^) and (g^, ^) form the 
complete set of equations which have to be solved to obtain the 2PI effective 
action at NLO in the expansion. From r[0, G] all correlation functions 
can then be found by derivatives with respect to the fields as functions of the 
known (p and G. We stress that the 2PI-1/N expansion is done on the level 
of the effective action. There are no further approximations involved on the 
level of the evolution equations. 

We note that Eq. ( |5.7| ) contains a double integration over the time contour 
C which can be inconvenient for numerical purposes. It turns out that it is 
possible to disentangle the nested integrations by exploiting the function K^. 
Convoluting the functions B and D and using the definitions for B and 
in Eqs. (E^) and (|5.2|), one obtains 



D(a;,2/) = i^5c{x -y) + ^Ka{y,x)(j)aix) 



j^Gabix,z)Gab{x,z)'D{z,y). (5.8) 



We observe that the nested integrals have disappeared in Eq. ( |^ ) without 
any problems. In Appendix B we work out more details for the equations 
preserving the nested-integral structure. It is also convenient to rewrite Eq. 
( |5.2| ) for Kq such that B does not appear. By convoluting B and Kq, one 
obtains 

A A /" 

^a{x,y) = —(t)b{x)Gba{x,y) -i— J Gbc{x,z)Gbcix,z)Ka{z,y). (5.9) 

As a result, we find that B and B~^ are eliminated completely from the 
coupled equations. We also see that the gap equations for D and are local 
in one time variable (here y), which is useful for numerical implementation. 

The form of the equation of motion for the inverse propagator, Eq. ( p.6| ), 
is suitable for a boundary value problem and can be used, in particular, to 
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discuss the propagator in thermal equihbrium by specifying the contour C 
to the Matsubara contour along the imaginary-time axis. However, to deal 
with nonequilibrium time evolution, i.e., an initial value problem, the form 
already given in Eq. (|2.13| ) is more useful. At NLO we find that the partial 
differential equation for the propagator takes the form 

+ + ^ + Gcc{x, x)] ^ Gab{x, y) = iSabSc{x - y) 

+(/)a(x)Kfe(x,?/) - i B{x,z)Gac{x,z)Gcb{z,y). (5.10) 



To summarize, Eqs. ( [5.1| , [5. 101 ) together with Eqs. ( |5.8| , p.9[ ) form the complete 
set of equations which have to be solved. They are completely equivalent to 
our first set containing Eq. ( WM for the inverse propagator. 



6 The auxiliary-field method 



In this section we show that the equations derived in the previous sections 
can also be obtained in the auxiliary-field formulation. In particular, we 
discuss the difference between the 2PI-1/N expansion and the "bare- vertex 
approximation" introduced in Ref. [pl|] . 

Following Refs. [0, ^ we rewrite the action by introducing an auxiliary 
field X as 



s[v,x] = - 



1 2N SAT 2 1 



(6.1) 



Integrating out x yields the original action ( p.l| ). From the Heisenberg 
equations of motion we see that the auxiliary field represents the composite 
operator x{.x) = \/ {6N)ipa{x)ipa{x). The following one- and two-point 
functions can be written down: 



and 



(paix) = {iPa{x)), X{x) = {X{x)), (6.2) 



Gab{x,y) = {Tc(Pa{x)ipbiy)) - {Vaix)){ipbiy)), 

Ka(a;,?/) = {Tcx{x)va{y)) - {x{x)){Va{y)) = 'Ka{y,x), (6.3) 
'D{x,y) = {Tcx{x)x{y)) - {x{x)){x{y))- 
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We note that additional "propagators" appear due to the introduction of the 
auxihary field. Since x is not a dynamical degree of freedom, only Gab has a 
physical meaning of a propagator. The role of the other one- and two-point 
functions is to implicitly perform infinite resummations, which were carried 
out explicitly in the direct method in the previous sections. 

Following the quantum fields are combined in an extended field with 
+ 1 components, 

^^ = {^:A^ (6-4) 



where i = l,...,iV-|-l. One may now formulate the 2PI effective action 
for this quantum field theory by coupling sources to the field ^i{x) and the 
bilocal field ^i{x)^j{y) ||2T|]. We would like to point out that in the presence 
of sources the original quantum theory and the one with the auxiliary 
field are potentially different. For instance, in the second formalism one 
may differentiate with respect to the bilocal source and obtain correlation 
functions of x{x)fa{y)- This possibility is not present in the original theory. 
However, as we will see below, on the level of the equations of motion in the 
2PI-1/N expansion at NLO, the two approaches yield identical results. 
The inverse of the classical propagator is 

where $j = {(pa,x) and S'[$] = S'[0, x]. It has the following components: 

S'S[<f>,x] 



S(l)a{x)6(j)b{y) 

S'S[(j>,x] ^ SN 
5xix)6xiy) A 

s's[cp,x] 



6x{x)6(j)a{y) 



- [D^ + + x{x)] SabScix - y), 

6c{x-y), (6.6) 
(j)aix)5cix -y). 



These operators are symmetric. Similarly, the matrix containing the two- 
point functions is defined as 

Gab -Kfj \ , _s 

K, D J' ^^-^^ 
where K.a{x,y) = Ka(?/, x) Hence, this matrix is also symmetric. 
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The 2PI effective action can now be written down and reads 

r[$, g] = + ^Tr iTiQ-^ + '-Trg^^g + V^iG] + const. (6.8) 

Here r2 is given by all two-particle irreducible graphs made with lines 
representing the "propagators" Gab, K^, K^, and D, and the vertex 
— \x{.^)Va{.x)ipa{x). In the auxiliary-field formalism, T2 does not depend 
on the expectation value $j. 

The equations of motion for the field expectation values follow by 
variation of F with respect to 0a and x- We find 

- [n^ + + x(x)] = ^ [K„(x, x) + K„(x, x)] = K„(x, x), (6.9) 

and 

X{x) = ^[<l>\^)+Gcc{x,x)\. (6.10) 

The equation for the two-point function follows by variation with respect to 
which gives 

e«' = eo:i-2i^. (6.11) 

By convoluting this equation from the right with Q and decomposing the 
self-energy as ^ 



- U n j . (6.12) 

one obtains the following set of coupled equations: 

- [D^ + + x{x)\ Gab{x, y) = 0a(a;)Kft(x, y) + i5abSc{x - y) 

+i J tac{x,z)Gcb{z,y) + Ea{x,z)Kh{z,y)^, (6.13) 

—Ka{x,y) = 4)b{x)Gba{x,y) 

+1 / [Eb{x, z)Gba{z, y) + Ii{x, z)\^a{z, y)] , (6.14) 

J z 

3N 

— D(x, y) = (f)a{x)Ka{x, y) + i5c{x - y) 



^Notc that in the auxiliary-field formalism Yiab does not receive a local LO contribution. 
It differs therefore from the self energy in the direct method, obtained by varying 
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Figure 6: NLO contribution in the 2PI-1/N expansion in the auxihary- field 
formahsm. The full line denotes the scalar propagator Gab, the dashed line 
the auxiliary-field propagator D. 

+t I pa{x,z)Ka{z,y) + Tl{x,z)Ti{z,y)\ , (6.15) 

J z 

- + + x{x)\ Ka(a;, y) = 0a(a;)D(x, y) 

+i tab{x,z)Ki,{z,y) + Ea{x,z)B{z,y) . (6.16) 

J z 

We note that Eq. ( |6.16| ) for is not an independent equation since 
Ka(a;, ?/) = K.a{y,x). Therefore, Eq. ( p.l6| ) is not needed in practice. 

To find explicitly at which order in the 2PI-1/N expansion specific 
diagrams contribute, we note that in the auxiliary-field formalism the possible 
irreducible 0(A^)-singlets are of the form 

tr(G"^), D, tr(KG'"K), (6.17) 

with n>l. From Eq. (|6rT5| ) it follows that D ~ 1/A^, and from Eq. ( |6?T^ 
we find that 

tr(KG'"K) ~ -^tii<P<pG"+^) (6.18) 

is proportional to 1/A^ as well. Using this scaling behaviour it is 
straightforward to give the diagrams that contribute at NLO and NNLO 
in the 2PI-1/N expansion in the auxiliary- field formalism. 

We find that the NLO contribution to r2 consists of one graph only, 

rr°[6^] = ^/ Gabix,y)Gabix,y)B{x,y), (6.19) 

shown in Fig. ^ From this expression the self-energies defined above follow: 
S^,L°(a;,y) = -Gab{x,y)-D{x,y), 

nNLO(x,i/) = -^Gab{x,y)Gab{x,y), (6.20) 
Sr(a:,l/) = 0. 
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a) 



b) 



c) 



Figure 7: NNLO contribution in the 2PI-1/N expansion in the auxihary- 
field formahsm. The full-dashed lines denote the mixed propagators K^, K^. 



Inserting these expressions in Eqs. (|6.13| )- (|6.15| ), we immediately recover our 
final result at NLO, Eqs. |5l^ , [5l^ , |5.1(]| ), obtained by the direct method 
in Sec. 

Only three diagrams contribute in the auxiliary-field formulation at 
NNLO. They are shown in Fig. |^. We note that diagrams with the mixed 
propagator (resulting in a nonvanishing Sq and Sq) appear only at NNLO. 



In Ref. 1 21 1 the first NNLO diagram, 



^NNLOa 



^ Jxy 



Kaix, y)Gab{x, Z/)Kfc(x, y), 



(6.21) 



was combined with the NLO diagram of Fig. ^ in the so-called bare- 
vertex approximation (BVA). We conclude that in the presence of a 
field expectation value the BVA approximation is not consistent with the 
2PI-1/N approximation discussed here. For vanishing 0, the 2PI-1/N 
approximation at NLO and the BVA ansatz are identical.^ 



7 Evolution equations for the spectral and statistical 
functions at NLO 

In order to describe nonequilibrium dynamics we will now specify the contour 
C to the standard Schwinger-Keldysh contour along the real-time axis ||24|] .p| 
The two-point function can be decomposed as 

GUx,y) = G>,{x,y)ec{x'-y') + G<,{x,y)ec{y'-x'), (7.1) 



^Therefore, we use the same notation for the functions K, D and 11 in Sects. |^ to ^ 
^For similar approximation schemes see also Ref. 

®We use a Gaussian initial density matrix. Non-gaussian initial density matrices are 
also possible, see e.g. Refs. 0]. 
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where G^^{x,y) = G^i^*{x,y) are complex functions. For the real scalar field 
theory it is convenient to express the evolution equations in terms of two 
independent real-valued two-point functions, which can be associated to the 
expectation values of the commutator and the anti-commutator of two fields 
i, I, i. We define 

Fab{x,y) = ^(G>,{x,y) + G<{x,y))=Re[G>,{x,y)], (7.2) 

Pab{x,y) = i(G>f,{x,y) - G<f,{x,y)^ = -2lm[G>f,{x,y)]. (7.3) 

Here F is the statistical propagator and p denotes the spectral function, with 
the properties F*f^{x,y) = Fab{x,y) = F^aiy.x) and p*b(x,?/) = pab{,x,y) = 
-pba{y,x). 

In order to proceed it is convenient to separate the singular part of D 
(see Eq. (|5^)) and write 

J:){x,y) = ^[i5cix-y) + -b{x,y)^, (7.4) 

with 

D(x, y) = Ka{y, x)0„(x) - ^n(x, v) + ^ J n(x, z)B{z, y). (7.5) 

For the functions Ka(x, y) (see Eq. (|5.9|)) and D(x, y) we define the statistical 
and spectral components as 

K{x,y) = ^(K>ix,y) + K<ix,y))=Re[K>ix,y)l (7.6) 
K^(x,i/) = i{K>{x,y)-K<ix,y))=-2lm[K>ix,y)], (7.7) 

and the same for I)F{x,y) and I)p{x,y). 

Now we have all the necessary definitions and relations to express the 
time evolution equations for the field expectation value and Green's function 
along the Schwinger-Keldysh contour as real and causal equations. The time 
evolution equation for the field reads 

- [D, + m^ + x{x)] M^) = (x, x) (7.8) 

with 

xix) = -^[4>\x) + FUx,x)]. (7.9) 
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The functions and satisfy the equations 

dzUp(x,z)K^(z,y) 

r dzUj.(x,z)KP^(z,y), (7.10) 
Jo 



X 



3iV 



^^ix,y) ^ :^Mx)pba{x,y) + — j ^ dzTlp{x,z)K.P^{z,y), 



where we employ the notation 



/ dz = dz^ I dz , 
Jo Jo J 



(7.11) 



and 



^F{x,y) = 



Fab{x,y)Fabix,y) - -pab{x,y)pabix,y) 



^p{x,y) = -Fab{x,y)pab{x,y). 



(7.12) 



A 



The statistical propagator obeys 

- [n^ + + x{x)] Fab{x, y) = -^Fac{x, x)Fcb{x, y) + 0a(a;)Kf (x, y) 

+ dztP^{x,z)Fcb{z,y) - dzt^^{x,z)pcb{z,y), (7.13) 
Jo Jo 

and the spectral function 

- [D^ + + x{x)] Pab{x, y) = :^Fac{x, x)pcb{x, y) + 0a(x)K^(a;, y) 



+ / dztp^^{x,z)pcb{z,y). 



Here we use the notation 



Kbi^^y) = 



A 



x_ 

3N 



(7.14) 



Fab{x, y)'DF{x, y) - -^pabix, y)t>p{x, y) , (7.15) 
Pab{x,y)t>F{x,y) + Fab{x,y)'bp{x,y) , (7.16) 
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with 

^3iV j ^^^p^^^^^^P^^^y^ ^ ^ J dzUF{x,z)'Dp{z,y), 

XX 

Dp(x,2/) = -KP^iy,x)(f)aix) - —np{x,y) + — I dzUp^x, z)'bp{z,y). 

(7.17) 

In absence of a field-expectation value {(f)a = 0) we find that Kf = = 
and the equations above reduce to those treated in 0, Q for diagonal two- 
point functions. 

In order to completely determine the time evolution, Eqs. ( |7.8| ), ( [7.13| ) and 
( 7.14 ) have to be implemented with initial conditions taken at x^ = y^ = 0. 
For the field (f)a{x) one may choose nonvanishing values (pai^^ = 0,x), but 
vanishing "velocities" dxO(j)a{x)\xO=o = 0. The initial values for the spectral 
function are completely fixed by the equal-time properties 

Pab{x,y)\xO=yO = 0, dxopab{x,y)\xO=yO = 6ab5'^{y^-y)- (7.18) 

Nontrivial information about the initial density matrix is contained in 
(derivatives of) the statistical two-point function at initial time 

Fab{x,y)\xO=yO=Q, dxoFab{x,y)\xO=yO=Q, dxodyoFah{x,y%o=yO=o. (7.19) 

Specification of these three functions is necessary and sufficient to solve the 
equations of motion. 



8 Weak coupling expansion 

In order to simplify the interrelated set of nonperturbative NLO equations 
of motion given in the previous Section, we discuss here a truncated version 
of the 2PI-1/N approximation for the case that the coupling is weak. It 
amounts to expanding the effective action to second order in explicit factors 
of the coupling constant A. Since Fg*"* is proportional to A it is preserved 
completely. For the NLO contribution F^'"^ in Eq. ( [4.3|) we find 

no] ^ ^^n,x.x) + .(A)^£n(x.,)n,x.,) 
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Figure 8: Diagrams contributing in the 2PI-1/N expansion when an 
additional weak-couphng expansion to second order in A is performed. 

-2z(^A^ j^U{x,y)Mx)Gab{x,y)Uy), (8-1) 

where we use again the notation Il{x,y) = —^Gab{x,y)Gab{x,y). The cor- 
responding diagrams are presented in Fig. 

The weak-couphng expansion affects the equations of motion through the 
auxihary variables and D, 

^aix,y) ^ ^(l)bix)Gba{x,y) + i (^-^^ j^U{x, z)(f)b{z)Gba{z,y), 

D(x,i/) - ^[Mx)Gabix,y)My)-^ix,y)]- (8.2) 

In the evolution equations (|7.13| ) and (|7.14|) for the statistical and spectral 
function the Eq. ( [7.10[ ) is replaced by 



A f X \ ^ 

Kf (x, y) ~ —(f)b{x)Ff,a{x, y)+ i^j / Up{x, z)(f)b{z)Fha{z, y) 

' A V r° 

J dzUF{x,z)(j)b{z)pba{z,y), (8.3) 



- -^<Pb{x)pba{x,y) + j j ^ dzY[f,{x, 



Z)^b{z)pba{z,y), 



^Here we consider the simple case of a field expectation value for which {X/6N)(j)'^ is 
small compared to the characteristic mass scale. We stress that for a consistent weak 
coupling expansion it is important to note that a nonzero minimum at (jp — 0q of the 



classical potential in (2.1) scales as 0q — —{QN/\)m?. Therefore, in a situation with 
spontaneously broken symmetry it is not sufficient to count only the powers of A coming 
from the vertices, as exemplified in this section. A consistent 0{\) scheme would have to 
take into account the first and the third graph of Fig. ^, while at 0{\^) the three-loop 
graphs of Figs. || and || would have to be taken into account. 
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and Eq. ( |7.17] ) simplifies considerably to 

Di.(x,y) ^ :^[(f)aix)Fabix,y)(j)biy) -IlF{x,y)], 

Dp(x, y) ^ ^ [(l)a{x)pab{x, y)(j)b{y) - np(a;, y)] . (8.4) 

9 Summary 

We have derived the 2PI effective action r[0, G] for the 0{N) model using 
the 2PI-1/N expansion to next-to-leading order. A detailed discussion of 
the classification of diagrams was presented. The equations of motion for 
the field expectation value (p and the two-point function G were calculated 
without further approximations. We showed the equivalence of the direct 
calculation with the auxiliary-field formulation. 

A detailed, but separate investigation would be necessary in order to 
discuss the question of the nonperturbative renormalizability of r[0, G] and 
the evolution equations derived above within the approximations considered 
in this paper. In principle, this problem may be treated following methods 
outlined in Ref. However, concerning the applications we have in mind 
and which are listed in the introduction, we emphasize that the physics of 
these problems is dominated by soft excitations and requires a finite cutoff. 
Therefore, from the practical point of view the important next step is to 
solve these equations along the lines of Refs. [|, |^ using a straightforward 
lattice discretization. 
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Appendix A 

In this section we discuss possible further approximations consistent with an 
expansion in powers of and comment on some aspects of Goldstone's 
theorem. We note that within the NLO approximation the replacement 

trG^ ^ (trG)VA^ (A.l) 
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in the expressions for the diagrams discussed in Sec. ^ is correct up to 
higher order terms. This corresponds to the replacement Gab{x,y)Gab{x,y) 

[Gaaix,y)]yN in the functions B(G) of Eq. Q and 1{G) of Eq. 
One observes that the resulting expressions can no longer be represented by 
the diagrams of Sec. |[ To verify this replacement up to NNLO corrections 
we note that (for given space-time coordinates) G can be diagonalized by 
virtue of 0{N) rotations. In particular, G is diagonal up to subleading 
corrections. This can be seen explicitly from the LO solution of Eq. ( p.9|) 
for the propagator, G^^^"* ~ Sab, which follows from the fact that the LO 
diagrams depend only on the invariants trG and (j)"^ (cf. Sec. ^). Since the 
invariant trG^ does not appear in LO diagrams, the replacement ( [A.l| ) is 
correct within the NLO approximation. We stress here that ( [A.l|) has not 
been used in the derivation of any equation presented in this paper. 

A similar argument cannot be applied to the invariant tr(0</)G) = 
4'aGab4'b- To see this, it is sufficient to restrict our attention to constant 
field configurations. In this case, it follows from 0{N) symmetry that the 
most general form of the propagator is 

gT\^) = GL{<P')Ptb + GT{<P')Plb. (A.2) 

where = (f>a(pb/<t>'^ and = 6ab — Pab are respectively the longitudinal 
and transverse projectors with respect to the field direction. Using this 
decomposition, we first check that the replacement ( [A.l|) is valid at NLO, in 



agreement with the above general discussion. Indeed the difference 



trG^ - = (Gl - GtY [l--]^N' (A.3) 

is subleading (recall that trC^ ~ trG ~ A^). However, 

tr(#G) - = <P\Gl -Gt){i-1-A-N (A.4) 



demonstrates that ii{(l)(j)G) cannot be replaced by (p'^tr G/N up to higher 
order corrections. 

In the remainder of this Appendix we want to show that Goldstone's 
theorem is fulfilled at any order in the 2PI-1/N expansion. Following Sec. ^ 
the 2P/ effective action can be written as a function of the 0{N) invariants 
(B only, 

r[0, G]=T [0^ tr(G"), tr(00G'P)] . (A.5) 
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In the case of spontaneous symmetry breaking one has a constant 7^ 
and the propagator can be parametrized as in Eq. ( |A.2|) . The standard IP I 
effective action riP7[(^] is obtained by evaluating the 2PI effective action at 
the stationary value (|2.9| ) for G 0, and the mass matrix A4ab can then be 
obtained from 

(52r[0,G(^*^*)(0)] 



M 



ab 



If rip7-[0] is calculated from ( |A.5D and ( [A .21 ) one observes that indeed the 
IPI effective action depends only on one invariant, 0^. The form of the mass 
matrix Aiab can now be inferred straightforwardly from rip/[0]. To obtain 
the effective potential ?7(0^/2), we write 



IP/ 



Ml 



=const 



(A.7) 



where Qd+i is the d+1 dimensional euclidean volume. The expectation value 



of the field is given by the solution of the stationarity equation (|2.7|) which 
becomes 

— — = 0a =0, (A.8) 

0(Pa 

where U' = dU / d{(jP /2) and similarly for higher derivatives. The mass 
matrix reads 



92f/(0V2) 



ab 



'^aO<Pb 



SabU' + (Pa(pbU" = iU' + 0't/")Pab + U' P^. (A.9) 

= 0) one finds that all modes have equal mass 
squared = U'6ab- In the broken phase, with 0^ 7^ 0, Eq. ( [A.8| ) implies 
that the mass of the transverse modes ~ U' vanishes identically in agreement 
with Goldstone's theorem. For a similar discussion, see Ref. ||2H1. Truncations 



In the symmetric phase 



of the 2PI effective action may not show manifestly the presence of massless 
transverse modes if one considers the solution (^(^^^t) of Eq. ( ^.91 ) instead of 
the second variation of r[0, G^^^^^^ for constant fields. For an early discussion 
of this point see Ref. as well as the comments in Ref. 



Appendix B 

In this Appendix we present the equations for the statistical and spectral 
functions preserving the nested integral structure and keeping the "chain of 
bubbles" I(a;, y; G) as the basic quantity. 
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All local contributions can be combined in an effective mass parameter 

M^,{X) = [W? + X{x)]5ab + ^ [<Pa{x)(Pb{x) + Gah{x, x)] , (B.l) 



and Eq. (|2.13| ) can be written as 



- \pa:5ac + MI^{x)\ Gcb{x, y) = iSab^cix ~ + ^acix, z)Gcb{z, v) , (B.2) 

with the "nonlocal" self-energy at NLO (we suppress the G dependence) 

T^abix, y) = {I(x, y) [(f)a{x)(f)b{y) + Gab{x, y)] + P(x, y)Gabix, y)} . 

(B.3) 

Here we defined 

P{x,y) = B-\x,u)A{u,v)B-\v,y), (B.4) 

= -(f)aix)Gabix,y)(l)biy). (B.5) 
Eq. ( [B.2| ) results in the standard time evolution equations for F and p ^ 

pJac + Mlix)] F^{x, y) = - dz z)F,,{z, y) 

Jo 



dz^acix,z)pcbiz,y), (B.6) 

PX" 

pJac + Ml{x)] pcb{x, y) = - dz S^,(a;, z)p,hiz, y), (B.7) 



with 



M^,{x) = [m^ + xix)]6ab + ^ [0a(x)0b(x) + Fab{x, x)] , (B.8) 
and the nonlocal self-energies 
^Ibi^^y) = -^{'^Fix,y)[(f)aix)My) + Fab{x,y)]-^Ip{x,y)pabix,y) 

+PF{x,y)Fab{x,y) - ^Pp{x,y)pab{x,y)y (B.9) 
Kbi^^y) = -^\M^^y)[^aix)(l)biy) + Fabix,y)] + lFix,y)pab{x,y) 



+Pp{x,y)Fab{x,y) + Pir(x, y)pafe(a;, (B.IO) 
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The functions Ip and Ip satisfy 

A 



A A f 

lF{x,y) = -—Upix^y) + — J dzIp{x,z)UFiz,y) 



3iV J, 



dz lF{x,z)Tlp{z,y), 



A A I 

Ip(a;,l/) = -^np(x,y) + — dzlp{x,z)Yip{z,y), (B.ll) 



and the nested integrals are 

A 



^F{x,y) 



AF{x,y) 

/•X" 

- / dz [Ap{x,z)lF{z,y) + lp{x,z)AF{z,y)] 
Jo 

+ / dz [AF{x,z)lp{z,y) + lF{x,z)Ap{z,y)] 
Jo 

- dz dvlp{x,z)AF{z,v)Ip{v,y) 
Jo Jo 

+ dz dvlp{x,z)Ap{z,v)lF{v,y) 
Jo Jo 



+ J dz j ^ dvlF{x,z)Ap{z,v)lp{v,y)Y (B.12) 



and 



Pp(a;,2/) 



A 



3iV 



Ap{x,y) 



dz [Ap(x, z)lp{z, y) + z)Ap{z, y)] 
+ dz dvlp{x,z)Ap{z,v)lp{v,y) 

JyO JyO 



(B.13) 



with AF{x,y) = -(j)a{x)Fab{x,y)(j)b{y) and Ap{x,y) = -(j)a{x)pab{x,y)(l)b{y)- 
The RHS in Eq. ( [T^ ) for the field expectation value reads 

^a{x,x) = :^Fab{x,x)(j)b{x) 
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3A^ 



dy \^p{x,y)Fab{x,y) + lF{x,y)pab{x,y)]<pb{y)- (B.14) 



Note that the nested time integrals in Eqs. ( |B.12| , |B.13|) have been ehminated 
in the equations discussed in Sects. ||, ^ by a convenient choice of auxihary 
variables. 
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